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Abstract. Lyman-break galaxies (LBGs) contain a non-negligible amount of dust. Flkkeuchi et al ] tOOSei) (T03) 
constructed a model of the infrared spectral energy distribution (SED) for very young galaxies by taking into 
account the dust size distribution in the early stage of galaxy evolution, which can be different from that of 
present-day evolved galaxies. We applied T03 model to LBGs and constructed their expected SED. In order to 
examine the grain size distribution of dust, we calculated the SEDs based on two distinct type of the distribution 
models: a single-sized distribution and a power-law distribution with a slope of dN/da oc a~^'^. We found that the 
single-sized and power-law dust size distributions yield a very similar detectability of LBGs at the submillimetre 
(submm). We also found that galaxies with a power-law dust distribution have much less flux at mid-infrared 
(MIR) than the other. By making use of this fact we can explore the dust grain size distribution in high-redshift 
galaxies through (observer-frame) FIR observations. Then, we applied the model to a gravitationally lensed LBG 
MS 1512— cB58 (cB58), a unique probe of the dust emission from LBGs. Observations by SCUBA suggest that 
the dust is hot in this galaxy. Our model framework well reproduced the hot dust temperature under a natural 
assumption for this galaxy. We also examined the detectability of LBGs at submm wavelengths in an eight-hour 
deep survey by ALMA. The LBG population with an age ^ 10** yr and a SFR ^ lOMQyr"^ can be detected in 
such a survey. By integrating over their redshifted SEDs with the observed luminosity functions, we obtained the 
contribution of LBGs to the cosmic infrared background radiation (CIRB). Although they have non-negligible 
amount of dust, their contribution was found to be small, especially in the FIR ~ 200 nm. Thus, we need a 
strongly obscured population of galaxies which contains a large amount of star formation, at some epoch in the 
history of the universe. 
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1. Introduction 

Star formation and metal enrichment in an early stage of 
galaxies are of great importance to understand the history 
of the universe. With the aid of a variety of new observa- 
tional techniques and large facilities, studies on such young 
galaxies are being pushed to higher and higher redshifts. 
The surveys that have most efficiently secured large popu- 
lation of galaxies at a redshift z ~ 3 are those that exploit 
the Lyman break to id entify high-z candidates: Lyman- 
break galaxies (LBGs) llSteidel et al.lll999tl2003l) . 

Even in LBGs, there is clear evide nce that they contain 
non-negligible amount of dust (e.g.. ISawicki fc ^^1^3119981 
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Adelberger fc Steidell lilOfit ICalzettil lioOlt IShaplev et all 
200lh . Dust grains absorb stellar light and re-emit it 
in the far infrared (FIR), hence it is crucial to eval- 
uate the intensity and spectrum of FIR emission from 
galaxies for understanding their star formatio n proper- 
ties (e.g.. lBuat et al.l2002tlHirashita et alJ2003l and refer- 
ences therein). Despite its importance, there still remains 
a large uncertainty in our understanding of the p roperties 
of du st emission from LBGs (for a review, see ICalzettil 

l200lh . 

lOuchi et alJ l|l999^ have estimated the submm flux 
from LBGs and, for the first time, found empirically that 
the dust temperature should be high in order to e xplain 
their non-detection by SCUBA. IChapman et alJ l)2000l) 
also pointed out that the SCUBA submm flux of LBGs 
were much weaker than that expected from their UV spec- 
tral slope index. They discussed that it is hard to recon- 
cile unless dust temperature T is higher than that of local 
starbursts (T k, 70 K). For the analysis of SCUBA flux 
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of a highly len sed LBG, MS 1512-cB58 (hereafter cB58), 
ISawickil 1I2OOII) (SOI) performed a thorough survey of the 
parameter space of dust temperature and emissivity index, 
and concluded that the dust temperature and/or emissiv- 
ity index in cB58 is substantially higher than those seen 
in local galaxies. 

Another important issue to be addressed is the contri- 
bution of dust emissi on from LBGs to the cosmic infrared 
background (CIRBI llGispert et al.ll2000HHauser fc Dwekl 



120011 iTakeuchi et all l2001at iTotani fc Takeuchil |2002|) 
The star formation rate (SFR) density of LBGs is now 
believed to increase t oward redshifts z ^ 3, if we 'cor- 
rect' dust extinction l|Steidel et al.l[l99^ . Then, the re- 
emission from dust in LBGs is an interesting target to 
examine. Some resent observations suggest that their con- 
tribution to the total CIRB is not very large, at most 
20 % ( Webb et a l. 2003). If so, there must be other con- 
tributors to the CIRB, which should bear an intense star 
formation comparable or even larger than that in LBGs. 
Hence, theoretical prediction of the dust emission plays 
an important role in the interpretation of the cosmic SFR 
and the CIRB. 

In order to address these issues and to make a con- 
sistent picture of the dust emission from LBGs, we must 
consider their star formation history, chemical evolution 
and dust mass evolution, radiation process with nonequi- 
librium temperatu re fluctuation, and d ust composition 
properly. Recently, iHirashita et al.l l)2002l) (H02) have de- 
veloped a model for the evolution of dust content in 
very young galaxies (age t ^ 10^ yr). Type II supernovae 
(SNe II) are the dominant source for the production of 
dust grains in young star-forming galaxies ( Dwek fc Scald 
H02 have modeled the evolution of FIR luminosity 
and dust temperature in such a young starburst on the 
basis of SNe II grain formation model of lTodini fc Ferraral 
1I2OOII) (TFOl). 

ITakeuchi etlll l)2003a(l (T03) subsequently con- 
structed a model of infrared (IR) spectral energy distri- 
bution (SED) of galaxies starting from the H02 model. 
It has been believed that the size of dust grains formed 
in SNe II cannot be as large as 0.1-1 /xm (TFOl), and 
host galaxies are too young for grains to grow in the in- 
terstellar space. T03, for the first time, properly consid- 
ered the dust size distribution peculiar to the very early 
stage of galaxy evolution, and construct a model of the IR 
SED of very young galaxies. T03 model successfully repro- 
duced the peculiar SED of a local metal-deficient galaxy 
SBS 0335-052 (Z = 1/41 Zq), which we regarded as an 
analogue of ge nuine young ga l axies. However, in clear con- 
trast, recently iNozawa et al.l ()2003 lh has presented a new 
model of dust grain formation in SNe II, which produce 
larger grains (^ 0.1-1 /im) than the previously believed 
conjecture. Their dust has a broken power-law grain size 
distribution. Based on the known physical condition in 
LBGs, we discuss the possible geometrical configuration 
of dust, and show that we can test the dust size distribu- 
tion at high redshift by an observation of LBGs in the IR 



and submm. Such observations will be soon possible by, 
e.g.,SPICA, HersciieJ, SMA and ALMA. 

The layout of this paper is as follows: In Section |2| 
we briefly describe our model framework. We present our 
basic results for the SEDs of LBGs in Sectional Related 
discussions are given in Section^ First we focus on a grav- 
itationally magnified LBG MS 1512— cB58, a unique probe 
of the internal physics of LBGs. Then we consider some 
cosmologically important problems. Section is devoted 
to our summary and conclusions. 



2. SED Model for Extremely Young Galaxies 

2.1. SED construction 

Since H02 treats the evolution of dust content in a galaxy 
younger than 10* yr, the only contributor to the total dust 
mass (Mdust) in a young (< 10* yr) galaxy is the supply 
from SNe II. The rate of SNe II is given by 



lit) 



^{t — im) <t>{in) dm , 



8M0 



(1) 



where ■0(t) is the SFR at age t (we define i = at the 
beginning of the star formation), 0(m) is the initial mass 
function (IMF), is the lifetime of a star whose mass 
is TO, and it is assumed that only stars with to > 8 Mq 
produce SNe II. H02 assume a constant SFR, -0 = V'O: 
and a Salpeter IMF. Then the rate of increase of Mdust is 
written as Mdust = Wdust7 (wdust is the typical dust mass 
produced by one SN II). H02 adopted TOdust = 0.4 Mq. 

In general we must take into account not only the pro- 
duction but also destruction of dust grains. However, the 
timescale of dust destruction is much longer than that we 
consider in this work, hence the effect of the destruction 
appears to be negligible. Detailed evaluation can be found 
in H02 [their eq. (2)]. 

Dust grains are roughly divided into silicate and car- 
bonaceous ones. TFOl have presented that the sizes of 
silicate and carbonaceous grains formed in SN II ejecta 
to be about 10 A and 300 A, respectively. Classical stud- 
ies claimed that the dust grains originating from SNe II 
cannot be as large as those formed in the atmosphere 
of evolve d AGB stars or grown in the diffuse interstel- 
lar space ijKozasa fc Hasegawalll987t TFOl). In this case, 
the discrete and small grain sizes make the appearance of 
the IR SED of young galaxies drastically different from 
that of aged n o rmal galaxies. On the other hand, recently 
iNozawa et all l|2003l) have proposed a drastically differ- 
ent picture of dust size distribution from SNe II. They 
showed that, under some condition, dust grains can grow 
large even within the expansion timescale of SNe ejecta. 
Consequently, their size distribution of grains have a bro- 
ken power-law shape. For comparison, we also consider 
this case by using t he simple description of Galactic dust 
llMa,this et a]lll977^ as 



dN 
da 
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for both grain species, i.e., silicate or carbon dust. Here 
we note that we do not include polycyclic aromatic hydro- 
carbons (PAHs) in our dust grain model, because they are 
rarely found in metal-poor systems and/or in intense ra- 
diatio n environments fe.g.. lMaddenll200ot iGalliano et al.l 
I200.1I) . 

It is well accepted that very small grains are stochasti- 
cally heated, that is they cannot estab lish thermal equilib- 
rium with the ambient ra diation field ( Draine fc Andersonl 
1198,'jlnraine fc Lill200lh : the heat capacity of very small 
grains is too small to maintain its temperature until the 
next photon impinges. As a result, the temperature of a 
grain varies violently in time, and we can define a distri- 
bution of the instantaneous temperature of a grain, which 
will be introduced in eq. ©• For the heat capacity of a 
grain, C (T) , we apply the multidimensional Debye model 
according to iDraineTLl ||2003). 

Then we calculated the temperature distribution of 
dust as a function of size a b y Mon te Carlo simulations, 
following IPraine fc Andersonl <| 19851) . The UV radiation 
field strength is determined by the history of OB star lu- 
minosity of II02 and the size of the star forming region 
rsF- For the geometry, we assume that the OB stars are 
concentrated in the centre of the system, and dust is as- 
sumed to be distributed as a shell. We regard the radius 
of the sphere, r as the size of the spherical star forming re- 
gion, rsF- The flux is expressed by J{t) = LoB(i)/(47rrgp), 
hence the UV energy density, u, is obtained by u = J{t)/c. 
The spectral UV energy density ux is calculated from the 
assumed spectrum of the UV radiation field. 

The rate at which a grain absorbs a photon with energy 
E ^ E + dE is expressed as^ 



dEdt 



(3) 



where Qabs is the absorptio n efficiency of a dust g rain. We 
used the values proposed bv lPraine fc Leel l(l984l) for Qabs 
of silicate and carbon grains. 

The heating is represented as follows: 



he 
T 



C{T')dT' , 



(4) 



where T is the peak temperature achieved by a grain hit 
by a photon with energy hc/\, and Tq is the grain tem- 
perature just before absorption. On the other hand, the 
grain cools through radiation as 



dT _ Stt 
H ~ ~ aC{T) Jo 



?abs(a,A) Bx{T)dX, 



(5) 



where Bx{T) is the blackbody function represented as a 
function of wavelength. 

From Equations ©,(01), and 10, we obtain the sample 
path of a grain temperature as a function of time. 



^ We note that there are two typographic errors in 
Equations (13) and (14) of T03. The equations we show here 
are the correct expressions. 



The total mass of each grain component is given by 
TFOl. The mass ratio we adopt here is Msh : Mc = 
0.56 : 0.44 (II02). With this value and material density 
of each species {p^» = 3.50gcm~'^ and pc = 2.26gcm~^: 
iDraine fc Le3ll984|) . we obtain the normalization of the 
dust size distribution dNi/da[, 



Mdnstfi — 



Aira'^p, dN, 
3 da' 



da'i 



(6) 



where subscript i denotes the species of dust, sil or C, 
and fi is the mass fraction of dust of species i. The total 
number of grains, Ni, is expressed as 



da'. 



da'i 



(7) 



For the TFOl dust distribution, is fixed for each dust 
component to a specific value a^, i.e., dNi/da[ — Ni5{a[ — 
a.i). In this case, with Equations lO and lO, the normal- 
ization reduces to 



3 A/dust /j 
47rafpi 



(8) 



We drop the prime in the following expressions. 

Total IR emission spectrum from a galaxy is then cal- 
culated by superposing each of the continuum from silicate 
and carbonaceous grains. 



W=E'^/ / 47ra?gi,,(A)S,(r) 



dN,dPi{u,a,) 



daj 



dT 



dTdtti (9) 



where dPi{u,ai)/dT is the temperature probability den- 
sity function, i.e., the fraction of time for which a grain 
with species i in a UV radiation field with energy density u 
stays in a temperature range [T, T+dT] . The starburst age 
is incorporated with Ni (i.e., dNi/dui) and dPi{u, ai)/dT 
through Afdust(i) and u\{t). Here we put a constraint that 
the amount of absorbed light is equal to the FIR lumi- 
nosity. Lastly, if the dust opacity is very large, then self- 
absorption occurs, and even MIR radiation from dust is 
absorbed by the dust itself. We treat the extinction as 
given by a shell model. Wc crudely approximate the dust 
opacity and extinction to the first order as 



7'dust(A) ~ 7ra^Qabs,i(A) 

i 

i 

/(A) = /o(A)e-^--(^) . 



3iV,; 



47r [(rsF + ArsF)3 - r|p] 

47rr2 ' 
*^'SF 



ArsF 
(10) 

(11) 



The absorbed light is re-emitted at longer wavelengths, 
mainly in the submm, and consequently, the SED is de- 
formed by the self-absorption. Final SED is obtained via 
this absorption-re-emission process. For further technical 
details, see T03. 
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2.2. Input parameters for LBGs 

Here we consider the input physical parameters in the 
model calculation for LBGs. We adopted the dust grain 
sizes ac — 200 A and = 6A for our canonical model. 
They are the same as those for SBS 0335—052 and I Zw 18 
used in T03, because LBGs are thought to be young, and 
hence we expect similar dust properties to them. 

The SFR of LBGs spreads over the range of SFR ~ 
1-300 M(r,yr-^ with a median value of SFR ~ 20 Mqjt-^ 
(e.g., lErb et al1l2003(l . We assumed a constant SFR up 
to the age shorter than 10^ yr, which may be a good ap- 
proximation (e.s.. ISaker et al1l2004l) . Within the young 
age of i ^ 10® yr, neither SNe I nor RGB/AGB stars con- 
tribute to the dust production. Other dust growth mech- 
anism in the in terstellar medi um cannot work effectively, 
either (see, e.g., IWhittetll 99?l. Thus, the basic framework 
of T03 model is vahd for LBGs. 

The most important information to calculate the IR 
SED is the effective size of the star forming region, but it 
is the most uncertain quantity at the same time. It is still 
beyond the possibility of the present facilities to measure 
the size directly from observations. The mean half-light 
radius of LB Gs is estimated to be ~ 1.6 kpc from HST 
observations llErb et al.ll200l . hence it may be a general 
value for LBG populations. The similarity between the 
rest UV and optical morphologies suggests that the dust 
in LBGs may not preferentially obscure particular regions 
in these galaxies l|Dickinsonl boOOt iGalzettil l200lh . Then 
we can safely use the galaxy radius vq as the radius of a 
star-forming region, rsF, i.e., rsF = 2 kpc. We will revisit 
this issue with more physical considerations in Section 0] 



3. Results 

Based on the above-mentioned settings, we calculated the 
IR SEDs of LBGs with SFR = 10, 30, and lOOMgyr"! as 
representative cases. We traced the evolution of the SEDs 
from t — 10^ yr to 10® yr. 

Figure^shows the evolution of the SED for the single- 
sized dust size distribution. The dust opacity is found 
to be not very large, and consequently, the silicate fea- 
tures appeared in emission. Only in the case of SFR = 
lOOMQyr^^, these features becomes a significant absorp- 
tion after 3 x 10^ yr. Most of the LBGs have a lower SFR 
than SFR — lOOA/oyr^^, hence basically we can expect 
the silicate emission band feature in their MIR spectra. 
Along the same line, since the self-absorption is not sig- 
nificant for LBG spectra, we also expect a strong MIR 
continuum emission from LBGs if the single-sized dust 
distribution takes place. 

For comparison, we next show the evolution of the SED 
for the power-law dust size distribution in Figure[21 In this 
case, an obvious difference is that the MIR radiation is an 
order of magnitude weaker than that for the single-sized 
dust results. This is because the power-law size distribu- 
tion has much smaller amount of small grains (a ^ 100 A) 
than the single-sized distribution. The extinction proper- 



ties are similar to that of single-sized dust case, despite 
the drastic difference of the SEDs. Thus, if the geometri- 
cal configuration of dust is similar among LBGs, the MIR 
observation can be a strong tool to explore the size distri- 
bution of dust grains. 



4. Discussions 

4.1. MS 1512- cB58 

It is still difficult to compare a model prediction 
for young galaxies to observational data of LBGs. 
Among the LBGs, a typical LBG cB58 is a fortu- 
nate exceptio n: it has been observed in various wave- 
lengths (e.g., [EUingson et alJlTflfla iRechtold et abllTflflTt 
iNak anishi et al.i il997[ iTeolitz et alj l2000t Baker et alJ 
200ljr 

because of a st rong ma gnifi catio n by gravitational 
lensing (factor 22-40: ISeitz et alJll99^ . We adopt a com- 
monly used value of 30 for magnification. The half-light 
radius of cB58 is estimated to be rg ^ 1 . 4 /1"^ k pc from a 
detailed image reconstruction ISeitz et a,1.l ill 9981) ^. 

From spectroscopic observations of nebular lines, the 
S FR of cB58 is est imated to be ~ 10-20h-'^ Mgyr"! 
IPettini et all feoool) . The age of major star forma- 
tion in cB58 was estimated to be ^sf — 35 Myr 
(jMattcucci & Pioino 2002). An older age es timate of 
tsF — 140 Myr is given by l|Baker et alJ l2004j) . and we 
adopt both of these values for our calculations. This is 
broa dly consistent with other ob servational suggestions 
(e.g.. lPettini et al.l2000t f Ellingson et alll99^ . and signif- 
icantly younger than the median value of those of general 
LBGs (- 350 Myr: Shaolcv ct al. 2001). 

The optical exti nction of cB58 is e stimated to be 
E(B - V ) - 0.27 JPettini et al.ll200d lEningson et"ai] 
Il996; Tc plitz et aLll2000|) . Emis sion from dus t in cB58 
is measured at two wavelengths: ISawickil ()200l|) observed 
this galaxy at 850 ^m by SCUBA on the JCMT, and no 
signal was detected above 3a level of 3.9 mJy- Ano ther 
observation by SCUBA l|van der Werf et al i bOQ l l) de- 

iBakerei 



tected a 850 /xm flux to be 4.2 ± 0.9 mJy. iB aker et alJ 
([200 1) reported a detection at 1.2 mm by MAMBO (Max- 
Planck Millimeter Bolometer) array, and the flux was 
1.06 ±0.35 mJy (4.4cr). 

The submm flux from the foreground cD galaxy (6" 
away from cB 58) may contaminate the observed data 
iSawickil I200I. After careful examination of the new 
submm mapping data, it has been concluded that the con- 
tamination i s ~ 5 % at 850 /zm and ~ 22 % at 1.2 mm, 
respectively (jBaker et alJ l2004V though there still remain 
some other uncertainties (e.g.. Baker ct al. 2004 discusses 
a possible contribution of emission lines or background 
sources. We may need an interferometric observation to 
solve this problem). In addition, the uncertainties in the 
lensing correction and in the UV based SFR estimation 
may also be as large as that of the contamination correc- 



^ Throughout this paper, we adopt a cosmology with 
(ft, no, Ao) = (0.7, 0.3, 0.7), where h = Hq/IOO. 
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1 10 100 1000 1 10 100 1000 1 10 100 1000 

Rest wavelength X [fj.m] Rest wavelength A [fJ.m] Rest wavelength K [/xm] 



Fig. 1. The evolution of the SED for the single-sized dust size distribution. From left to right, the star formation rate 
(SFR) is 10, 30, and lOOA/oyr^^. The dotted, dashed, dot-dashed, long-dashed, and solid lines represent that the 
age of the major star formation in a galaxy is 1.0 x 10^ yr, 3.0 x 10^yr,1.0 x 10* yr, 3.0 x 10* yr, 3.0 x 10^ yr, and 
1.0 x 10^ yr, respectively. 




1 10 100 1000 1 10 100 1000 1 10 100 1000 

Rest wavelength A [fJ-m] Rest wavelength A [/^ni] Rest wavelength A [^m] 



Fig. 2. The evolution of the SED for the power-law dust size distribution. Comparing with Figure ^ we see a clear 
difference in the MIR. 



tion. We should keep these issues in mind in the following 
discussions. 

4.1.1. IR SED 

We use a conrmonl y accepted value of 24 M0yr~^ as a SFR 
ijBaker et al.ll2004|) and assume the SFR to be constant in 
time. For the age of the major star formation in cB58, 
as we mentioned above, we adopt 35 Myr and 140 Myr. 
Because of its young age, dust is predominantly produced 
by SNe II and the framework of T03 model is again valid 
for cB58. The ha lf-light radius of t he reconstructed image 
of cB58 (~ 2knc lSeitz etHlllflQi^ is very similar to that 
of general LBGs l.Bkpc lErb et alJ l2003'). 

We show the model SED of cB58 in Figure El 
Left panel presents the SED with a single-sized dust, 
and Right panel shows the SED with a power-law dis- 



tribution (Equation [^J, respectively. Syr nbols re p resent 
the measured f lux densities repor ted by ISawickil l|200lh 
(upper limits), 'Baker et al.| l)200lf) (open triangles) and 
Ivan do r Werf et al.. (,,2001) (open squares). There are two 
open symbols at each wavelength: the upper open symbols 
are the measured values, and the lower ones represent the 
flux densities corrected for the possible cont amination by 
the p oint source at the lens cluster center ijBaker et alJ 
l2004j) . It is still not clear whether the point source is the 
cD galaxy or not (Baker et al. 2004). 

The apparent peak of the SED is located at a wave- 
length shorter than 100 /im at the restframe of cB58. 
Actually, the dust temperature is higher than that of 
the local (z ~ 0) sample of normal galaxies observed by 
SCUB A, but comparable to the most intense dusty star- 
bursts (jPunne fc E ales"200l'). It is mainly because of an 
intense UV radiation field of cB58 and low dust opac- 
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Fig. 3. Model spectral energy distribution (SED) of a Icnsed Lyman-break galaxy MS 1512— cB58 (cB58). Left panel 
shows the observed SED calculated based on the single-sized dust size distribution with silicate grain radius asii — 6A 
and carbon grain radius ac = 200 A. The observed SED is magnified by a factor 30. Right pan el shows the SE P with 
a power-law grain size distribution. Symbols represent the measured flux densities reported bv lSawickill)200l|) (upper 
limits), leaker ct al. (2001) (open triangles) and van der Werf et al. (2001) (open squares). In both panels, Dotted lines 
represent the SED with the age of tgp — 35 Myr, and dashed lines of ^sf — 140 Myr, respectively. 



ity. The mean dust temperature of the power-law dust 
model is lower than that of the single-sized dust model, 
but the superposition makes the submm part of its fi- 
nal SED similar to that of a hot dust emission. We find 
that both models predict slightly lower flux densities than 
observed. Considering the uncertainties mentioned above, 
we conclude that the model is roughly consistent with ob- 
servations within a few factors. Some authors claim that 
the submm flux d ensities of cB58 is surprisingly low (e.g., 
iBaker et al.ll200-3) . but in our framework, they are natu- 
rally explained by our model framework (the model pre- 
dicts even smaller fluxes). The power-law model predicts 
a slightly larger flux at submm wavelengths, but it is diffi- 
cult to distinguish the two models only by submm obser- 
vations. 

The prominent feature in the MIR is the silicate band 
of 9.7 ^m. For the single-sized dust model, the NIR-MIR 
part of the SED is dominated by the radiation from small 
silicate grains, and FIR-submm part by larger carbona- 
ceous grains. In contrast, relatively lower abundance of 
small grains in the power-law model makes the MIR radi- 
ation considerably weaker. 

4.1.2. Extinction and size of the star-forming regions 

Besides IR emission, dust extinction is a fundamental 
quantity to specify the dust properties. We relate the ex- 



tinction in magnitude {A\) with Tdust(A) (Equation I10|l 

by 

Ax - 1.086Tdust(A) . (12) 

When we adopt the tsf = 2 kpc, the extinction at X-band 
is Apc = 0.08 mag. A gain we adopt the extinction curve of 
ICardelli et al.l ()l989l) . we expect Ay = 0.8 mag from our 
model. The measured color excess is E{B — V) = 0.27, 
hence by ad opting Ry = 3.1 as a representative value in a 
diffuse ISM llKnige]ll2003|> . we have Ay = RyE{B -V) ~ 
0.84 mag, which is in p erfect agreement with our model 
prediction. IViih et all ()2003l ) performed a sophisticated 
analysis of the attenuation in LBGs by using a radia- 
tive transfer model. They showed the attenuation is in 
the range of 5. 7-18.5 at 1600 A, an d is well represented 
by Calzetti law (jCalzetti et al.lll994(l . Our extinction esti- 
mate is also consistent with their result. They also showed 
that the shell geometry is appropriate for the dust in 
LBGs, which is in accordance with our assumption. Thus, 
we see that the extinction in the visible stellar system of 
cB58 is consistent with its dust amount and spatial scale 
of rsF — 2 kpc. 

Then, can a heavily hidden star formation exist in this 
galaxy? Actually, SB S 0335—052 hosts a completely ex- 
tingu ished starburst (jHunt et alJl200ll: IPlante k. Sauvagj 
I2DO2!), and its SED is well understood by T03 model. To 
consider this problem, we have calculated the SEDs with 
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Wavelength X [fJ-m] 



Fig. 4. Effect of the size of the star-forming regions on 
the model SED of cB58. Sohd Une is the canonical model 
(same as the solid line in the left panel of Fig.|3J). Dotted, 
dashed, and dot-dashed lines indicate the SED calculated 
with a radius of the star-forming region, rsp = 300 pc, 
120 pc, and 30 pc, respectively. Note that the SED is pre- 
sented in the unit of vL^, [Lq] here. 

various smaller star- forming region sizes. This means that 
the bulk of the dust in cB58 could be confined in a small 
part of its volume. They are shown in Figure 0] Dotted, 
dashed, and dot-dashed lines are those with rsF = 300 pc, 
120 pc, and 30 pc, respectively. Smaller sizes make the 
opacity in the region significantly larger, and the con- 
version from UV to FIR becomes much more efficient. 
Considering the observational upper limits, the allowed 
smallest size is rsp = 120 pc, with a possible range of 
50-200 pc if we take into account the uncertainties in the 
estimates of the SFR and lensing magnification. 

It has been pointed out that the galactic wind can 
transport the dust int o larger radiu s (T03) . Actually, pre- 
cise spectroscopy by IPettini et al.l l)2003|) clearly shows 
the outflowing motion of gas with a velocity of Wwind — 
255 km s^^ in cB58. They consider the wind and extinc- 
tion, and suggest that the superwind shell size is compa- 
rable to the galaxy size. Hence the wind crossing tiniescale 
is icross ~ 10^ yr. If dust closely couples with the gas mo- 
tion, it is difficult to keep dust grains confined in a re- 
gion much smaller than the galaxy size through its age of 
t = 35-140 Myr. Strong radiative pres sure rtiay al so ex- 
pand the dusty region, as suggested by (llnoue!'200^. Here 
we should comment the effect of the patchiness of dust, 
which may be realistic for LBGs. In the FIR, since the 
opacity of dust is small (t < 1), the FIR emission escapes 
from a galaxy without significant absorption. Then the 



flux is controlled simply by energy conservation, and the 
patchiness does not affect the result s ignificantly (see, e.g., 
IWitt et aljEooiillCordou et a,].ll2000h . 

Together with its high SFR and SN rate, which are 
considered to provide a large amount of radiative and ki- 
netic energy to gas-dust system, we conclude that a com- 
pact dusty star formation may not be plausible for cB58. 
Hence, the dust configuration of cB58 may be similar to 
the present model. This suggests that we can safely dis- 
criminate the effect of dust size distribution on the SED 
from that of the configuration and geometry of dust in the 
galaxy. Again, this conclusion can hardly be affected by 
the dumpiness of dust, just the same reason mentioned 
above. 

4.1.3. Follow-up observations of cB58 

Our single-sized dust model provides an SED with bright 
continuum in the rest NIR-MIR regime, hence we ex- 
pect that cB58 will be detected by MIR-FIR observations 
(A ~ 20-100 //m). At present, Spitzer Space Telescope^ 
MIPS and ASTRO-F^ IRC will be suitable instruments 
for such observations. For example, both instruments have 
the same level of sensitivity around ~ 20 /zm, and the fiux 
of 0.06-0.08 mJy wiU be detected. For FIR, both can de- 
tect 0.8-1 mJy at A = 60-70 /zm. Therefore, the MIR 
continuum can be detected by the future MIR-FIR obser- 
vations if the single-sized dust takes place 

On the other hand, if the power-law distribution takes 
place, the flux densities will be an order of magnitude 
smaller than those of our canonical model, because the 
number of the stochastically heated grains is reduced. As 
we mentioned, the uncertainty caused by dust geometry 
is not very large for a LBG because the configuration be- 
comes more or less shell-like in these galaxies just as a 
scale-up of SBS 0335-052. Therefore, the MIR observa- 
tion will be a strong test for the size distribution. 

We should note, however, that even by Spitzer. con- 
fusion limit can be severe in FIR (e.g. , Ilshii et aP 120021: 
iDole et al.ir2003HTakeuchi k IshiilbOO^ . Larger space IR 
facilities such as Herschel Space Observatory^ or SPICA^, 
and IR interferometry missions are longed to improve the 
confusion limit. At longer wavelengths, of course ALMA 
will be useful. 

4.2. Observability of the dust emission from LBGs 

Considering the non- negligible extinction in LBGs, their 
energy budget radiated in the IR s hould be significant 
(■Adelbereer fc Steidcl 200QclCalzettii200 1). Stimulated by 
the consideration, the relation between LBGs and dusty 
starbursts, another representative star-forming galaxy 
population at z = 2-5, has gathered the researchers' at- 

^ URL: 'http : //s irtf . caltech. e du/l 

URL: http : //www . ir . isas . ac . jp/ASTRO-F /index-e . html| . 
^ URL: http://www.rssd.esa.int/herschel/ . 
^ URL: http : //www . ir . isas . ac . jp/SPICA/index . htiiil| . 



8 



T. T. Takeuchi & T. T. Ishii: Dust Emission Model of LBGs 



tention. In fact, however, most of t he LBGs have not been 
detected by SCUB A observations I Chapman et alJl200Ct 
Webb et al.l| 2003), exc ept for a few extrem e cases (e.g., 
Chapman et al.. .2002(1 . ICalzetti et al.l l)2000|) pointed out 
the possibility that the dust is hot in LBGs based on the 
observed SED of a low-metaUicity galaxy Tol 1924—416. 
Hot dust in LBGs was also suggested by SOI from his 
observation of a lensed LBG cB58. Since their consider- 
ation is empirical, and we need some explanation for the 
hot dust, if it exists. As seen in the above, our model 
naturally predict a high dust temperature. Here, we dis- 
cuss how LBGs look like in the IR/submm wavelengths. 
It is also useful to estimate the observational feasibility of 
LBGs by forthcoming submm facilities like ALMA. 

The observed flux density of a source at observed fre- 
quency, z/obs is obtained by 



(1 + z)L(^i+^)^^^^ _ (1 + 



(13) 



where di^{z) is the luminosity distance corresponding to 
a redshift z, and Vohs and i^cm are observed and emitted 
frequency, respectively. We show the observed IR/submm 
SEDs of LBGs at z = 2, 3, and 4 in Figures EUTUl 
Figures |S}{7| represent the SEDs for the single-sized dust 
size distribution, and Figures 151 1 101 for the power-law dis- 
tribution (Equation 13) . The thick black short horizontal 
lines indicate the S-cr detection limits for 8-hour observa- 
tion by ALMA. Here we assumed 64 antennas and three 
wavelength bands, 450, 850, and 1080 /im. We also show 
the 3-cr source confusion limit of Herschel at 75, 160, 
250, and 350 /xm bands by thick gray horizontal lines. 
These limits are bas e d on 'the p hotometric criterion' of 
I Lagache et al.lEo03l ). Ilshii et al.l |2002) also report simi- 
lar estimates for the confusion limits.^ 

From Figures ISl 1 101 we recognize that the detectability 
of these galaxies at submm bands are not significantly de- 
pendent on the redshift. This is because of the 'negative 
X-correction', well-known in the field of submm astron- 
omy. On the contrary, the galaxy age and SFR are more 
important for the detection of LBGs. If the age ^ 10* yr 
and SFR ^ 10 Mq yr, a LBG can be detected at a wide 
range of redshifts in the submm. 

In general, the longer the wavelength is, the easier the 
detection becomes for millimetre (mm) observations. In 
the submm and mm wavelengths, the detectability is not 
quite different between LBGs with single-sized dust and 
those with power-law dust. It is a natural consequence 
that the radiation in these wavelengths is dominated by 
relatively large dust grains that can establish a thermal 
equilibrium with the radiation field. 

The apparent peak of the dust emission lies in a rel- 
atively short wavelength of about 200-400 /im at z = 
2-4. Thus, the accumulated background radiation spec- 
trum from LBGs may have its peak in the FIR. The ex- 



Generally, in a sky region suitable for cosmological surveys, 
the cirru s confusion is les s severe than the galaxy confusion 
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(see, e.g.,|Doiee£a0|2Q03). 



Fig. 11. The contribution of LBGs to the cosmic in- 
frared background radiation (CIRB) (open symbols). 
Filled symbols an d hatched area are measured CIRB spec- 
tru m by COBE (iHauser et all Il998t iFixsen et all Il998l: 
Lagache et alJll999(l 



pected dust emission from LBGs is, however, too faint 
to be detected even by Hershel deep survey, except some 
rare exceptionally bright sources, e.g., lensed galaxies like 
MS 1512— cB58, or in the case of extremely high SFR (in 
order to reach the confusion limits, unrealistic integration 
time might be required). Hence, still more powerful in- 
struments at these wavelengths would be required. 

In the shorter wavelengths (MIR~FIR) , the size distri- 
bution affects the possibility of detection drastically, just 
as we have seen in the case of cB58 (|4.1(l . Surveys in the 
MIR by Spitzer and ASTRO-F will be desirable to detect 
these objects. These observations are important to com- 
pare the overall SEDs of LBGs directly with those of "local 
analogues", such as SBS 0335—052 and I Zw 18. Existence 
of such local galaxie s with hot dust is reported by re- 
cent st udies (see e.g.. lCalzetti et allbOOnt iTakeuchi et alJ 
l2003bl) . Similarity and difference between them provides 
us a unique clue to the physics of galaxy formation. 

4.3. Contribution of LBGs to the CIRB 

The implication of hot dust also offers a potentially im- 
portant cosmological insight into high-z galaxies: if dust in 
high-z object is as cool as usually assumed, then we cannot 
expect a high comoving SFR density because in a realis- 
tic cosmology the cosmic IR background ( CIRB) spectrum 
strong ly constrains the high-z IR emission llTakeuchi et alJ 
|2Q^3l- However, very hot dust in such high-z object 
reconciles the CIRB constraint with high SFR, and al- 
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Fig. 5. Prediction for the observed IR/submm SEDs of LBGs at z = 2, 3, and 4. The star formation rate (SFR) is 
10 Mq yr~^ , and the dust grains are single-sized. The dotted, dashed, dot-dashed, long-dashed, and solid lines represent, 
the same as Fig.Q that the age of the major star formation in a galaxy is l.Ox lO^yr, 3.0x 10^yr,1.0x 10® yr, 3.0 x 10® yr, 
3.0 X 10'' yr, and 1.0 x 10^ yr, respectively. Confusion limits of Herschel and detection limits of ALMA 8-hour survey 
are also shown by gray and black thick horizontal lines, respectively. 




Fig. 6. The same as Figure but the SFR is 30 Mq yr~^ 




Fig. 7. The same as Figure HI but the SFR is 100 Mq yr ^ . The silicate feature finally appears as absorption. 
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Fig. 8. Prediction for the observed IR/submm SEDs of LBGs at z ~ 2, 3, and 4. The star formation rate (SFR) is 
lOM0yr~^, and the dust grains are power-law. 



„ 100.0 

>^ 

Co 

10.0 







2=2.0 


SFR= 30.0Msyr-' 






// 

h 
//- .- 


A 
\ -\ 

\ \\ 

\" "A 

■■■ ' A 


f ' 


■■■■ \\ 


if: / 




' :' 


\ 

\ 

\ 




10 100 1000 

Observed wavelength A [/im] 



10 100 1000 

Observed wavelength \ [/^m] 



10 100 1000 

Observed wavelength X [/im] 



Fig. 9. The same as Figure |H1 but the SFR is 30 M© yr~^ 



lows vigo rous star formation hidde n by dust in the early 
Universe (|Totani fc Takeuchill2002|) . 

Actu ally, non-detection of LBGs in the subm m wave- 
lengths ()Chapman et al.ll2000l: IWebb et al.ll2003ll suggests 
the hot dust temperature. Considering the non-negligible 
extinction in these galaxies, their energy budget radi- 
ated in the IR is significant (lAdelberger fc Steidell Eooil 
ICalzettilEooTL However, their contribution to the submm 
radiation cannot be quite large because of the high dust 
temperature. In fact, the contribution to the extragalactic 
submm back ground radiation h as been suggested to be at 
most ~ 2n% (IWehb et a,lJl2nn3^ . 

Usi ng our model SEP and the luminosity function of 
LBGs tehaplev et al we calculate the integrated IR 

light from LBGs to the CIRB. In order to have a rough 
estimate, we correct the extinction to the observed lumi- 
nosity function, and converted it to the distribution of the 
SFR by converting the UV luminosity to the SFR accord- 
ing to our model. Here we crudely assumed that all the 
LBGs have the same star forming region size tsf ~ 2 kpc. 
We also assumed that the number of LBGs rapidly de- 



creased at z ^ 4.® By using the SFR distribution, we 
obtained the multiband FIR luminosity function of the 
LBGs. We finally integrated it over the FIR luminosity to 
get the CIRB intensity. 

Our model predicts 0.4nWm^^sr^^ at 200 /i, and 
0.2nWm"^sr-i at 850 /im, respectively. The difference of 
the dust size distribution does not affect the CIRB spec- 
trum almost at all. This result is shown in Figure [TTl The 
850 /im intensity is quite consis tent with the LBG contri- 
bution at 850 /im estimated bv lWebb et all l)2003|) . Since 
the CIRB intensity at 200 /tm is much smaller than the 
measured peak (vly = 15nWm~^sr~^ at A ~ 140 /im) 
iLagache et alJll99^ . it strongly supports the existence 
of heavily obscured dusty galaxy population. Therefore, 
LBGs are not the dominant population contributing to 
the CIRB. 

Straightforward conversion of the UV radiation from 
galaxies into the IR b y recent models can not reproduce 
the strong peak fe.g.. iBalland et allbOOSj) . In addition, 
the large energy budget radiated in the IR requires an ef- 



Actually, this assumption affect the results little. 
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fective energy release, which may related t o starburst phe- 
nomena, even if they are nearby galaxies llTakeuc hi et a.lJ 
l2001at iFranc eschini et alJ lioOlt iHauser fc DwelJ I2OOI'). 
Therefore, hidden starbursts are necessary to explain 
the strong evolution of submm source counts (e.g., 
iTakeuchi et al.ll2001allJ) . At present, redshifts of the star- 
burst population contributing to the CIRB are almost un- 
known. If the dusty starburst is a low-2 phe nomenon, there 
must be a strongly obscured era at z ~ 1 ijTakeuchi et alJ 
2001al) . In this case it is interesting to discuss a possi- 
ble relation between the activation of star formation and 
the peak of merge r rate in hi erarchical structure forma- 
tion scenarios fe.g. JLacev fc Cole-1993 : Kitavama fc S uto 
On the other hand, if the obscured galaxies reside 
in a high-z Universe, it implies that the high-z cosmic 
SFR densitv mav b e higher than the dust-corrected value 
l)Steidel et a,lJll99<l> because it would be a sum of the con- 
tributions both from LBGs and hidden starbursts. 

5. Summary and conclusion 

Lyman-break g alaxies (LBGs) contain a non-negligible 
amount of dust. FSkeuchi et al ] ^|2n03a^ (T03) constructed 
a model of the infrared spectral energy distribution (SED) 
for very young galaxies by taking into account the dust size 
distribution in the early stage of galaxy evolution, which 
can be different from that of present-day evolved galaxies. 

In this work, we applied T03 model to LBGs and con- 
structed their expected SED. In order to examine the grain 
size distribution of dust, we calculated the SEDs based 
on two distinct type of the dist ribution models: a single- 
sized distribution predicted bv iTodini fc Ferraral l|200l|) . 
and a power-law distribution with a slope of dN/da oc 
a~3-5^ which is oft en used to describe the Galactic dust 
l)Mathis et a]|ll977|) . 

We found that the single-sized and power-law dust size 
distributions yield a very similar detectability of LBGs at 
the submillimetre (submm). We also found that the dif- 
ference of the grain size distribution affects the SED dras- 
tically at mid-infrared (MIR) wavelengths. Galaxies with 



power-law dust distribution have much less flux at MIR 
than the other. Generally, strong outflow is observed in 
LBGs, and consequently their geometry of dust configura- 
tion might be relatively simple. It reduces the uncertainty 
of complex radiative transfer, and we can safely explore 
the dust grain size distribution high-redshift galaxies by 
(observer-frame) FIR observations. 

Then, we applied the model to a gravitationally lensed 
LBG MS 1512-CB58 (cB58), a unique probe of the dust 
emission from LBGs. Thanks to the large magnification 
factor, the dust emission has been detected in the submm 
by SCUBA observations. These observations suggest that 
the dust is hot in this galaxy. Our model framework well 
reproduced the hot dust temperature under a natural as- 
sumption for this galaxy. 

We also examined the detectability of LBGs at submm 
wavelengths in an eight-hour deep survey by ALMA. 
The LBG population with an age ^ 10^ yr and a SFR 
^ lOM0yr~^ can be detected in such a survey. The 
detectability of LBGs in the submm is not strongly af- 
fected by redshift, because of the well-known negative K- 
correction. Hence, the detected LBG sample will be un- 
doubtedly an ideal sample to study the early evolution of 
metal and dust content in galaxies. 

By integrating over their redshifted SEDs with the 
observed luminosity functions, we obtained the contribu- 
tion of LBGs to the cosmic infrared background radia- 
tion (CIRB). Although they have non- negligible amount of 
dust, their contribution was found to be small, especially 
in the FIR ^ 200 /im. Thus, we need a strongly obscured 
population of galaxies which contains a large amount of 
star formation, at some epoch in the history of the uni- 
verse. 
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